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The ampitude of the B — > py decay induced by the flavour-changing neutral currents contains 
the penguin contribution and the weak-annihilation contribution generated by the 4-quark operators 
in the effective Hamiltonian. The penguin contribution is known quite well. We analyze the weak- 
annihilation which is suppressed by the heavy-quark mass compared to the penguin contribution. 

In the factorization approximation, the weak annihilation amplitude is represented in terms of 
the leptonic decay constants and the meson-photon matrix elements of the weak currents. The 
latter contain the By, py transition form factors and contact terms determined by the equations of 
motion. We calculate the B'y and py form factors within the relativistic dispersion approach and 
f*"^ ■ obtain numerical estimates for the weak annihilation amplitude. 

o 

£ ■ I. INTRODUCTION 

' 

The investigation of rare semileptonic B decays induced by the flavour-changing neutral current transitions b — ► s 
and b — > d represents an important test of the Standard Model or its extentions. Rare decays are forbidden at the tree 
level in the Standard Model and occur through loop diagrams. Thus they provide the possibility to probe the structure 
of the electroweak sector at large mass scales from contributions of virtual particles in the loop diagrams. Interesting 
information about the structure of the theory is contained in the Wilson coefficents in the effective Hamiltonian which 
£f) \ describes the b — > s, d transition at low energies. These Wilson coefficients take different values in different theories 

■ with testable consequences in rare B decays. 

Among rare B decays the radiative decays b — > sy and b — > dy have the largest probabilites. The b — > sy transitions 
are CKM favoured and have larger branching ratios than the b — > dy transitions. The b — > sy transition has been 
observed by CLEO in the exclusive channel B — > K*y in 1993 and measured inclusively in 1995. The B — > py decay 
, will be extensively studied by BaBar and BELLE. 

The main uncertainty in the theoretical analysis of B decays is connected with long-distance QCD effects arising 
*-£h I from the presence of hadrons in initial and final states. In inclusive decays these effects are under better control, 
however from inclusive measurements it is more difficult to obtain precise results. 
O || The decay amplitude contains two different contributions: one arising from the electromagnetic penguin operator 
1} . and another from the 4-fermion operators in the effective Hamiltonian. One of the effects generated by the 4-fermion 
' operators is the weak annihilation (WA). Further details about short- and long-distance effects in the radiative decays 
can be found in recent publications [Q-Q and references therein. 

In the B — * K*y decay the weak annihilation is negligible compared to the penguin effect: it is suppressed by two 
powers of the small parameter A ~ 0.2 of the Cabbibo-Kabayashi-Maskawa (CKM) matrix. In B — > py both effects 

■ have the same order in A and must be taken into account. 
The penguin contribution has been analyzed in several ways and is known quite well. On the other hand, the 

WA in B — > py has been studied in less detail: the relevant form factors were analyzed within sum rules ^|J6) and 
perturbative QCD 0. However, some contributions to these form factors were neglected. These contributions may 
be relevant if precise measurements become available. 

The aim of this paper is to analyze the weak annihilation for the B~ — > p~y decay more closely. 
In the factorization approximation, the weak-annihilation amplitude can be represented as the product of meson 
leptonic decay constants and matrix elements of the weak current between meson and photon. The latter contain 
the meson-photon transition form factors and contact terms which are determined by the equations of motion. The 
photon can be emitted from the loop containing the b quark which is described by By transition form factors. It can 
also be emitted from the loop containing only light quarks described by the py transition form factors (Fig[l]). 

We calculate the By form factors within the relativistic dispersion approach which expresses these form factors in 
terms of the B meson wave function. We demonstrate that the form factors calculated by the dispersion approach 
behave in the limit m& — > oo in agreement with perturbative QCD. The B meson wave function was previously tested 
in the B — > light meson weak decays and is known quite well, allowing us to provide reliable numerical estimates for 
the By form factors. 
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FIG. 1. Diagrams describing the weak annihilation process for B — > p7 in the factorization approximation: (a) The photon 
is emitted from the loop containing the b quark, (b) The photon is emitted from the loop containing only light quarks. 

The p7 transition form factor is related to the divergence of the vector and axial-vector currents. In the case of the 
axial- vector current it is proportional to the light-quark masses if the classical equation of motion is applied. Because 
the quark momenta in the loop are high, these masses have to be identified with current quark masses. For this reason 
the corresponding p7 form factor was neglected in previous analyses [QJ^] . We find however that this argument is not 
correct: this form factor remains finite in the limit of vanishing light quark mass m — ► and behaves like ~ M p f p /M% 
which means the violation of the classical equation of motion. We present the result for the pj transition form factor 
but leave a detailed discussion of the anomaly appearing for the matrix element (p7|c? ! ,A„|0) for a special publication. 

Finally, we provide numerical estimates of the weak-annihilation amplitude taking into account the Bj, pj transition 
form factors, and the contact term contributions. 

In Section II the effective Hamiltonian for the b — > d transition and the general structure of the amplitude are 
presented. In Section III we discuss the photon emission from the B meson loop and obtain the £7 transition form 
factors within the dispersion approach. Section IV contains results for the pj transition form factors. In Section V 
the numerical estimates are given. The concluding Section summarises our results. 



II. THE EFFECTIVE HAMILTONIAN, THE AMPLITUDE AND THE DECAY RATE 

The amplitude of the weak radiative B — > p transition is given by the matrix element of the effective Hamiltonian 
for the b — ► d transition 

A(B pry) = ( 1 (q 1 )p(q 2 )\H cS (b -> d)\B(p)), (1) 

where p is the B momentum, q 2 is the p momentum, and qi is the photon momentum, p = qi + q 2 , q\ = 0, q\ = Mr, 
p 2 — M|. The effective weak Hamiltonian has the structure |^]: 

H oS (b -> d) = -| Z t C 7l (p)0 7l - ^j=t u (CiOu)Oi + C 2 ( P )0 2 ) , (2) 

where only operators relevant for the penguin and weak ahnnihilation effects are listed. Gf is the Fermi constant, 
£ g = V* d V q b, Ci's are the Wilson coefficients and 0,'s are the basis operators 

g _ 

077 = ttt d a <Ju.vrn b (ii) (1 + 75) b a F^, (3) 

01 = d a J v (l ~ l5)u a U/97„(l - 75)^) 

02 = rfa7^(l - 75)"/3 Upj v (l - 7 5 )fr Q , 

(4) 

with the following notation: e = \/ATra cnl , 7 15 = ij ^ 1 ^ 2 ^ 3 , cr M „ — i [7^,71/] /2, e 0123 = —1 and Sp (^y^^ ^) = 

The amplitude can be parametrized as follows 

cG 

A(B~ -> p-7) = —pr- [e 9lC . 92£ -Fpc + it*2t\ (.9^Wi - P^Qiu) F PV ] , (5) 

where Fpc and -Fpv are the parity-conserving and parity-violating invariant amplitudes, respectively. e 2 (ti) is the 
p-meson (photon) polarization vector. We use the short-hand notation e a bcd = f-ap^vO^b^ 'c^d" for any 4-vectors 
a, 6, c, d. 

For the decay rate one finds 

T(B- - p-j) = ^^M% (1 - Ml/Mlf (|F PC | 2 + |Fpv| 2 ) . (6) 
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A. The penguin amplitude 

The main contribution to the amplitude is given by the electromagnetic penguin operator C>7 7 : 

A pcng {B -> pi) = -^^^C 7 ^|Ti(0) (e gie * 92£2 « +i4 v H fl (9vnP<li ~ PnQiv)) > (7) 

where Xi is the form factor of the B — > p transition through the tensor current [p|-pl| . The corresponding contribution 
to the invariant amplitudes is therefore 



^ P P c nS = ^p P v S = -fcCV^lKO). (8) 



B. The weak annihilation amplitude 

The radiative B — > p 7 transition also receives contribution from the 4-fcrmion operators 0\ and O^. For the 
charged B~ — > p~ (92)7(91) transition the corresponding amplitude reads 

A WA (B- -» p-j) = --^|^„(p(<ft)7(gi)|d7»/(l " 7s)w • ii7„(l - l 5 )b\B(p)), (9) 

In what follows we suppress the lable WA. Neglecting the nonfactorizable soft-gluon exchanges, i.e. assuming vacuum 
saturation, the complicated matrix element in Eq. (^ is reduced to simpler quantities - the meson-photon matrix 
elements of the bilinear quark currents and the meson decay constants. The latter are defined as usual 

(p(q 2 )\d^u\0) = e^M p f p , f p > 0, 
{Q\vrf v jsb\B(p)) = ipufg, f B > 0. (10) 

It is convenient to isolate the parity-conserving contribution which emerges from the product of the two equal-parity 
currents, and the parity-violating contribution which emerges from the product of the two opposite-parity currents. 

1. The parity-violating amplitude 

The parity-violating amplitude has the form 

A PV (B pi) = ^ u a 1 {{pr/\d'y v u\0)(0\vrf v 'y B b\B) + (p\d^u\0)(iu 7ul5 b\B)} . (11) 

Here a\ is an effective Wilson coefficient, which we take as a± = C% + CijN c at the scale ~ 5 GeV. 

The first term is a contact term which can be calculated using the equations of motion for the quark fields [EJ. 
Setting m u = ma, we find 

(pid^u\0)(0\u lljl5 b\B) = ip v f B (rY\d7uu\0) = f B ((r/\du(d lv u)\0) = iee^e^ g^Mpfpfe- (12) 
The B — ► 7 amplitude in the second term of (O) can be parametrized as follows 



(l{<a)\uiv-1bb\B{p)) =iee*/ 



, ,2F A 2f B 



(13) 



It contains the form factor -Fa (92 = ^p)i ano - ^ ne contact term proportional to f B (the derivation of this relation is 
given in the Appendix). 

Summing the contributions of the photon emission from the i?-meson loop and the p-meson loop gives the amplitude 
^4pv which can be represented in the form Apy = t^ L T p with 9^ T), = as required by gauge invariance. Thus, the 
weak-annihilation contribution to the form factor Fpy for the B~ — > p~"f decay isn 



1 We note that for the B + — > p + 7 the term proportional to fs changes sign (see Appendix), and also F% = —F^ as can 
be obtained from charge conjugation of the amplitude Apy. For the B° — > p°7 decay the term ~ fs is absent. 
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F^ V A = Zuaif P M p 



2F A 



Ve 



M b M% - Mj 



(14) 



The two contact terms which are present in the amplitudes of the photon emission from the B meson loop and from 
the p-meson loop do not cancel each other (we disagree here with the claim of Ref. ||) but lead to a nonvanishing 
contribution proportional to fs- 



C. The parity-conserving amplitude 



This amplitude reads 



Gf, 



Apc (B -> pr,) = -^6.01 {{p\diuu\0){i\u 7u b\B) + ( 7 p\dw 5 u\0){0\uw 5 b\B)} . 
The B — > 7 amplitude from the first term in the brackets contains the form factor -Fy(g| = M^): 



(j(qi)\ui v b\B(p)) = ei 



M f 



(15) 



(16) 



The second term in ([15]) can be reduced to the divergence of the axial- vector current and contains another form factor, 
Gy\ namely, 



(01^75615) (7/91^75^10) = /s(7Pl^c?7^75"|0) = e f B G v e qi<i * q2f :* 
Thus, the weak annihilation contribution to -Fpc reads 

2-FV Jb Gy 



Fpc — ^uaiM p fp 



M F 



Mpfp 



(17) 



(18) 



Summing up, within the factorization approximation the weak annihilation amplitude can be expressed in terms of 
the three form factors Fa, Fy, and Gy. 



III. THE FORM FACTORS F A AND F v 



In this section we derive the formulas for the form factors Fa .v within the dispersion approach to the transition 
form factors. This approach has been formulated in detail in [|l 0|1 and applied to the weak decays of heavy mesons in 
p|. Recall that the form factors Fa,v describe the transition of the £?-meson to the photon with the momentum q\, 
qj — 0, induced by the axial- vector (vector) current with the momentum transfer q2, q\ = M^. For technical reasons, 
it is convenient to treat the form factor Fa(v) as describing the amplitude of the photon-induced transition of the 
B-meson into a bu axial- vector (vector) virtual particle with the corresponding factor l/(s — g|) in the dispersion 
integral. Then we can directly apply the equations obtained in J10| for the meson-meson transition form factors. 



A. The form factor Fa 



The form factor Fa is given by the diagrams of Fig |[ Fig ^a shows F^ , the contribution to the form factor of the 
process when the b quark interacts with the photon; Fig ||b describes the contribution of the process when the quark 
u interacts while b remains a spectator. One has 



F A = Fi« 



F 



(u) 



(19) 





iy 5 z — r^YvYs 



a. b. 

FIG. 2. Diagrams for the form factor F A - a) F A b) , b) F A u) . 
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The B meson is described by the vertex b(kb) «75u(£v) G(s) / \fW c , with G(s) = <Pb(s)(s — M B ). The S-meson wave 
function <f>s is normalized according to the relation [|10| 

oo 

^2 / ds(j) B {s) (s ~ {m b - m u ) ) = 1. (20) 

(m5+m u ) 2 

Here A (a, b, c) = (a — b — c) 2 — 46c is the triangle function. 

It is convenient to change the direction of the quark line in the loop diagram of Fig[^b. This is done by performing 
the charge conjugation of the matrix element and leads to a sign change for the 7^75 vertex. 

Now both diagrams in Fig |2] a,b are reduced to the diagram of Fig || which defines the form factor F^\mi, 1TI2): 
Setting mi = nib, m 2 = m u gives F^: — QbF^' (mb,m u ). Similarly, setting mi = m u , m 2 — nib gives F^\ 
F [ a ] = -Q u F i A 1) (m u ,m b ). 




^5 - k, m, XYs 



FIG. 3. The triangle diagram for '(7711,7712). The cut corresponds to calculating the imaginary part in the variable p . 

For the diagram of Fig || (quark 1 emits the photon, quark 2 is the spectator) the trace reads 

- Sp ij 5 (m 2 - ^2)7^75(^1 + &i)7n(wi + hi) = 4i(&i + k[) ^(m^ + m 2 k 1 ) v + Ai(g^q a - 3 MQ ^)(m 1 fc 2 + m 2 fci) Q . 

The spectral density of the form factor F^ (rn\, m 2 ) in the variable p 2 , p = k\ + fc 2 , is the coefficient of the structure 
g^u obtained after the integration of the trace over the quark phase space. Performing necessary calculations, we 
arrive at the following single dispersion integral 



Mb' 



-F 



(i) 



/Nr. 



At: 2 



(m b +m u ) 2 



ds 4>b{s) 
(«-M p 2 ) 



1 fX^(s,mlml) 2 



pqi 



2s 



ni\ log 
m 2 log 



s + m\ — m 2 + A x / 2 (s, m 2 , m|) 
s + m 2 — m 2 — A 1 / 2 (s, m 2 , m 2 ,) 

s + mf — m| + A 1 / 2 (s, m 2 , m|) 
s + m 2 — m 2 — A 1 / 2 (s, m 2 , m 2 ,) 



+ (m 2 — mi) 



A 1 / 2 (s,m 2 ,m 2 ) 



(21) 



For g 2 = one has pq x = (M 2 - M 2 )/2. 

Now, let us analyze the behaviour of the form factor in the limit mj —►00. To this end it is convenient to rewrite 
the spectral representation (|l]) in terms of the light-cone variables as follows (see JlJ] for details) 



2 (i), . \/N~ c f dx x dx 2 dk\ <f> B {s) , 2 -, 

— — 1\ (mi, m 2 ) = —5- / 5 0(1 — x\— x 2 ) 77^ 4 mnz; 2 + m 2 xi - (mi - m 2 )k ± pqi \ 

M B 47r z J x\x 2 s — M£ 



(22) 



Here Xi is the fraction of the B-meson light-cone momentum carried by the quark i, and s = m 2 / 'xi+m 2 ,/ 'x 2 +k 2 _/ 'x\x 2 . 



For the form factors F^ and FjP one obtains 



— M u) = -Q u - 
M B a 

2 f ( V - n 



N c f dxdkj_ 0s(s) 



4tt 2 
4tt 2 



x 2 Xfc s - M 2 
dxdk\ 4>b{s) 
%b^ u ^ M„ 



m u x b + m b x u + 



m b x u + m u x b + 



2(mfc — m u )k 



M% - Ml 
2(m u - m b )fc 2 
M% - M 2 



with 



mt m' 



(23) 
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Let us recall that the B-meson decay constant has the following representation in terms of the wave function |K 

\/W c f dxdk\ </>b(s) 



{m u x b +m b x u } . (24) 



47T 2 J X u Xb s — M? 



p 



Due to the wave function <f>g(s), the integral in the heavy quark limit is dominated by the region x u — A/mb, 
x b = 1 — A/mb, where A is a constant of order Mb — m b . This leads to the following expansion of the form factors in 
the l/m;, series 



M B a Am b 

J_ F W = Q A + (25) 
Mb mi 

Clearly, the dominant contribution in the heavy quark limit comes from the process when the light quark emits the 
photon, whereas the emission of the photon from the heavy quark gives only a l/m^ correction. The expressions ( pq) 

for the form factor agrees with the result of Ref. [Q , while we have found a different sign for ^ 



B. The form factor Fy 



The consideration of the form factor Fy is very similar to the form factor F^. Fy is determined by the two diagrams 
shown in Fig [|: Fig ||a gives Fy \ the contribution of the process when the b quark interacts with the photon; Fig [I|d 
describes the contribution of the process when the quark u interacts. One has 

F V = F$ ) +F V U) . (26) 




a. b. 

FIG. 4. Diagrams for the form factor Fy- a) Fy , b) F } 



It is again convenient to change the direction of the quark line in the loop diagram of Fig ^b describing F v by 
performing the charge conjugation of the matrix element. For the vector current 7„ in the vertex the sign does not 
change (in contrast to the 71/75 case considered above). 

Then both diagrams in Fig ^ a, b are reduced to the diagram of Fig || which gives the form factor Fy (mi, TO2): 
Setting mi = m&, m 2 = m u gives F^': F^' = QbFy (mb,m u ); Setting mi = m u , m 2 = m b gives F^, Fy" 1 = 
Q u Fy\m u ,m, b ). 




FIG. 5. The triangle diagram for Fy (mi, 7712). The cut corresponds to calculating the imaginary part in the variable p 2 . 



G 



The trace corresponding to the diagram of Fig 4 (1 - active quark, 2 - spectator) reads 

- Sp ij 5 (m 2 - k 2 )^ v {m 1 + fc'i)7M( TO i + = -4e yAIQ9l (mifc 2 + m 2 k 1 ) c 
The spectral representation for the form factor takes the form 



M B ' 



4tt 2 



(m b +r. 



ds(f) B (s) (, . A 1 / 2 (s,m?,m 2 ) / s + m\ - m\ + A 1 / 2 ^, m 2 , to 2 ) 

i^Mf) \ (m2 " Wl) s + mi l0g t S + TO 2 -TO 2 -AV 2(s , m 2 , m 2 ) 



To analyze the heavy quark limit to& — > oo we again represent the form factor in terms of the light-cone variables 



M B ' 



-F, 



(i) 



N c f dx 1 dx 2 dk 2 t </>b( s ) , , 
0(1 - xi - x 2 ) 775- (toix 2 + m 2 xi) . 



4tt 2 



x\x 2 



s-M 2 



(27) 



For i*v>" and K> the corresponding expressions read 



(u) 



v ^ /" dxdfc 2 ^ <j) B (s) 



x*x b s - M 2 



~T7 V ~ ~Qb- 



47T 2 J X 2 

v^/Vc /" dxdkj_ <j) B (s) 



An 2 J x\x u s - M 2 
By the same procedure as used for Fa, we obtain in the limit to?, — > 00 



{m& x„ + m„ x b } , 



Mr 



F< u > = -Q u 



^F {b} - -O h — 

s r V ~ 2 

M s v mi 



fb 



Ami, 



(28) 



The dominant contribution in the heavy quark limit again comes from the process when the light quark emits the 
photon. Now both form factors Fy^ and Fy^ in ( f28|) perfectly agree with the expansions obtained in 0). 

As seen from Eqs . ( |25| ) and (pq), one finds Fa — Fy in the heavy quark limit, in agreement with the large-energy 



effective theory 1 13 



IV. THE FORM FACTOR Gi 



The form factor Gy is determined by the divergence of the matrix element of the charged current between the 
vacuum and the p~j state, 



iPv{l{qi)p (<?2)|d7 I ,75tt|0) = eG v e qie * q2e * 
The corresponding diagrams are shown in Fig ra. 



(29) 





FIG. 6. Diagrams describing the amplitude ( , ~f(qi)p~ (q2)\djufsq\0), p — qi + 92, p 2 = M%. The diagram (a) is multiplied 
by the d-quark charge, and the the diagram (b) is multiplied by the u-quark charge. The cut corresponds to the calculation of 
the imaginary part in the variable gf- 
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If the classical equations of motion are applied, the form factor Gy isproportional to the light-quark masses. This 
is the reason why this form factor was neglected in previous analyses However, a proper calculation shows that 

this argument is not correct: in fact the classical equations of motions do not hold and the divergence contains the 
anomaly. 

The anomalous behavior of the divergence of the axial- vector current in the chiral limit is a well-known phenomenon 
discovered in the two-photon amplitude (77|S„g7„75<3'|0) (m). A very clear way to demonstrate the anomaly is to 
start with the matrix element of the axial vector current and to calculate the spectral representations for the relevant 
form factors. The anomaly is then obtained by performing the divergence at the final stage of the calculation p5| . A 
similar treatment applied to the matrix element (7P~ \d^ v ^^u\Q) leads to the form factor Gy which does not vanish in 
the limit m — > 0. We will present a detailed discussion of the anomaly in radiative B decays in a separate publication 
p6| . Here we only quote the final result: considering the spectral representation in the variable q 2 and introducing 
the appropriate p-meson radial wave function ip p (s), one finds the following expression for the form factor Gy 



G\ 



Ml r ds f p {s) (s 
4tt 2 J (s 



Ml 



M 2 ) 



(30) 



ipp(s) is normalized according to the relation (for massless quarks) 



oo 

^ J ds a |W(s)| 2 = 1. 



(31) 



Clearly, Gy is finite for m — which means a violation of the classical equations of motion for the axial- vector current. 
Eq. d30| ) describes the anomaly which takes place for the 7/5 final state as well as for the 77 one. There is however an 
important difference between the two cases: For the 77 final state the divergence remains finite in the limit m — > 
and p 2 — M% — > 00. For the p7 final state the divergence is finite for m — ► but decreases as 1/p 2 for p 2 — > 00. It is 
convenient to introduce the parameter k such that 



Gy = {Q u + Q d )n Mpfp 



Ml ' 



(32) 



with k staying finite for m = and Mb — > 00. 



V. NUMERICAL ESTIMATES 



Let us write once more the penguin and the weak annihilation contributions to the B — > p 7 amplitude: 

,pcn g _ F pc„ g = _ eta7 mfc Ti(0)j 

~2F A , 2/b 



PV 

WA 
PV 



PC 

iuaiMpfp 



M B Ml 



Ml 



-^PC 



Z-udxMpfp 



2F 



v 



JbG\ 



M b Mpfp 

The scaling behavior of the form factors in the limit Mb — * 00 reads 



Ti(0) ~ M B 3/2 [17], Fa,b - M\ 



-1/2 



such that 



F pong - M, 



-1/2 



F 



-WA 
PV,PC 



M 



G\ 



-3/2 



Ml 



(33) 



(34) 



(35) 



The terms proportional to Jb (/b ~ l/yMs) are l/M^-suppressed compared with the terms containing Fa,v- As 
we see below numerically this leads to a suppression by a factor of 4 — 5. 

We now proceed to numerical estimates for the B-meson decay. The scale-dependent Wilson coefficients Ci (p) and 
ai{jj) take the following values at the renormalization scale p ~ 5 GeV B: 



C x = 1.1, 



C 2 = -0.241, 



C 77 = -0.312, 



ax = d + C 2 /N c 



1.02. 



(36) 
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The penguin form factor was previously calculated within the dispersion approach with the result T x (0) = 

0.27 ±0.3 jllj]. Using the same parameters and the B meson wave function as in jll]] we obtain the form factors Fa,b 
shown in Table I. Our result for the form factor Fy is in good agreement with the estimates from other approaches. 
The form factor Fa agrees well with the constraints from perturbative QCD and turns out to be considerably larger 
than the corresponding sum rule estimate. 

The value of the Gy is very sensitive to the details of the p meson wave function ip p . The reason for that is the 
presence of the term (s — MT) in the integrand in Eq. ( |30| ) which changes sign in the integration region. Assuming 
tpp(s) ~ /J 2 /(/3 2 + s) 2 and setting (3 = 0.8 GeV, which gives a good description of the p meson radius, leads to 
k = —1.8. However, this value has a large uncertainty. Conservatively, we take |k| < 2.0 and use this result for further 
estimates. The form factor Gy does not contribute more than a few % to the full amplitude, but can sizeably correct 
the weak-annihilation part. 

Using the obtained form factors and the decay constants fg = 0.18 GeV, f p — 0.21 GeV we arrive at the following 
estimates (we show contributions of different terms in Eq. (p3|) separately): 



F po„ g = ^ 20 MeV; 

Fpy = £ u {-7.8 (cont. of Fa) + 2.2 (cont. of f B )} MeV = -5.6 f u MeV 

= f u {-6.0 (cont. of F v ) ± 0.7(cont. of G v )} MeV = -(6.0 ± 0.7) £ u MeV. 



(37) 



Taking into account errors in the form factors, for the ratios of the weak-annihilation to the penguin amplitudes 
we find 



parity — violating : 
parity — conserving : 



i^A /F Peng „ _ (0 28 ± . 025) 

F^c A /Frc S - "(0-3 ± 0.05) W6- 



(38) 



\£ u /£t \ ^ 0.4 in the Standard Model. 

Our results for F pcng and F WA agree with the sum rules § which reported the ratio F WA /F pcns ~ -0.3 
We would like to notice, however, that the anatomy of F WA in our analysis is different. For instance, we have found 
Fa considerably larger than the sum rule result. But after including the contact term ~ /g which was not taken into 
account in ||, we have come to Fpy^ close to the sum rule estimate. 

The calculated form factors and the ratios of the weak-annihilation to the penguin amplitudes is one of the necessary 
ingredients for the calculation of the branching ratios of the B — > pj decays, Isospin and CP Asymmetries. In addition 
to the above ratios, these quantities contain the phase induced by the strong interactions and the CP-violating phase 
of the CKM matrix (see for details). The corresponding analysis was done recently in || using the value 
F WA /FP cn s ~ -0.3 f tt /&. 



TABLE I. Results for the weak-annihilation form factors Fa, Fv and Gv for the B — * p 7 decay. The accuracy of our 
estimates is about 10%. The sum rule results are recalculated from |(| according to the relation Fa,v ~ —F^Mb / f p - ■ The 
results from j?j are recalculated according to Fa,v = -^/Jl for A = 0.5 GeV. 





Disp Approach 


SR p| 


pQCD § 


-F A 


0.120 


0.073 


>0.09 


-F v 


0.092 


0.091 


>0.09 


\Gv\ 


<0.004 











VI. CONCLUSION 



We have analyzed the weak annihilation for the radiative decay B — > pry in the factorization approximation. 

1. We have calculated the form factors Fa and Fy describing the photon emission from the B meson loop within the 
relativistic dispersion approach. We have performed the 1 /nib expansion of the form factors and demonstrated that 
the form factors of the dispersion approach exhibit a behaviour in agreement with the large-energy limit of QCD. 

2. We have analyzed the contribution to the weak annihilation amplitude from the diagram when the photon is 
emitted from the loop containing only light quarks. For the parity-conserving process this quantity is related to the 
divergence of the axial- vector current 

(7(91)^^(92)1^^75^10) = ee 91£ j 92e .(Q n + Q d ) n f p M p /M%. (39) 

with k staying finite in the chiral limit m — > 0. This result means the violation of the classical equations of motion 
and represents an anomaly which has the same origin as the anomaly of the matrix element (7(51)7(52)1^^10)- The 
value of k is found to be sensitive to subtle details of the p- meson wave function. Conservatively, we estimate |«| < 1. 

3. We have also included contact terms which were missed in some of the previous analyses. Numerical estimates for 
the weak annihilation contribution to the B~ — > p~j amplitude are given in Eq. d38|). 

It was noticed in || that the weak annihilation mechanism is crucial for the D — ► pry decays in which it dominates 
over the penguin mechanism. It is therefore very important to take into account all the contribitions to the weak 
annihilation amplitude listed above for a proper description of the rare radiative D decays. 
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VII. APPENDIX: TRANS VERSITY OF THE PARITY- VIOLATING AMPLITUDE. 



The parity-violating amplitude ( pi] ) is given by the sum of the two terms 

A PV (B~ - — - aPV 1 aPV 



PI 



A™+A 2 



(40) 



where A PV - {p~ (q 2 )\d lv u\0)( 1 (q 1 )\u ll/l5 b\B- (p)) and A PV = (p- {q 2 ) 1 {q 1 )\d llJ u\Q) (Q\u lvTo b\B- {p)) . 
1. We start with A pv . Let us write (j(qi)\uj u j 5 b\B- (p)) = ee^T^ with 

T*(p,q)=i [ dxe^(0\T(J^x),u lul5 b)\B-(p)), 



(41) 

where Jp(x) — |w7 M u— ^d^pd — ^bj^b is the electromagnetic quark current. The amplitude has the following Lorentz 
structure 

T »u = * (g^PH - P^qiu) F 1A {ql) +i (q\pp~pqi qi M ) q\ v F 2A (q\) + i (q\ Pp - pqi qip) PvF 3A (ql) + ~ ( 42 ) 

where C is the contact term. The contact term can be determined using the conservation of the electromagnetic 
current dp J M = 0, which leads to the relation 



iA (p, q) = - (0| [Q, <*y„7fc«] \B~ (p)) 



ifsPv = -ifBPv 



(43) 



for the B meson. This gives C = —Jb- Notice that for the B°-meson the contact term is absent. For the radiative 
decay q\ = so we find 



-4i V = ie fpMpe^e™ \ {g^pqi - Pp,qiv)Fi A {0) - PpPuj-. 



'21 



B 



iefpMpt^e™ < {g^pqi - Pp,qi„)F 1A (0) - p^qi 



M% - Ml 

Vb 
v Ml - Ml 



L B 



(44) 



2. Using the equation of motion for the quark fields (Q u — 2/3 e, Qd — —1/3 e) 

ilvd v q(x) = mq(x) - Q q A v 'fq{x), 
id v q(x)^ v = -mq(x) + Q q A u q(x)Y, 



one obtains for A P v 



A 2 =iPvfB{p 7|d7„u|0) = f B (p j\dv(dy v u)\0) 



3. For the sum wc find 



A\ v + A P v = ie fpfsMpt'/el" (g^pqi - Ppqi 



(Qd - Qu)e (p -f\A u d^u\0) = ee^e^g^yfpMpfB + 0{m u ,m d )- 

Vb 



F 1A (0) + 



M% - M* 



(45) 



(46) 



(47) 



and obtain the relation (|l4| ) after setting F\ A = 2F a /Mb- For the B° — > p°7 decay the term proportional to /a is 
absent. 
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